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[57] ABSTRACT 

A system for indicating ice thickness and rate of ice 
thickness growth on surfaces is disclosed. The region to 
be monitored for ice accretion is provided with a reso- 
nant surface waveguide which is mounted flush, below 
the surface being monitored. A controlled oscillator 
provides microwave energy via a feed point at a con- 
trollable frequency. A detector is coupled to the surface 
waveguide and is responsive to electrical energy. A 
measuring device indicates the frequency deviation of 
the controlled oscillator from a quiescent frequency. A 
control means is provided to control the frequency of 
oscillation of the controlled oscillator. In a first, open- 
loop embodiment, the control means is a shaft operated 
by an operator. In a second, closed-loop embodiment, 
the control means is a processor which effects auto- 
matic control. 

8 Claims, 14 Drawing Figures 
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MICROWAVE ICE ACCRETION METER 

GOVERNMENT CONTRACT 

The invention described herein was made in the per- 
formance of work under NASA Contract No. NAS3- 
22765 and is subject to the provisions of Section 305 of 
the National Aeronautics and Space Act of 1958 (72 
Stat. 435; 42 USC 2457). 

DESCRIPTION 

1. Technical Field 

The present invention relates to an apparatus for 
indicating ice thickness and rate of ice thickness growth 
on surfaces which may be subject to ice accretion as a 
result of weather conditions. The invention finds partic- 
ular utility in the avionics field for detecting and mea- 
suring ice thickness and rate of ice growth on air foil 
surfaces. 

2. Background Art 

The problems caused by icing, particularly in the 
avionics field, and some solutions to those problems 
heretofore provided by the prior art are summarized in 
my prior U.S. Pat. Nos. 4,054,255 entitled “Microwave 
Ice Detector” and 4,060,212 entitled “De-Icing Appara- 
tus and Method”. The research which led to the inven- 
tions described in the above-mentioned patents is more 
thoroughly discussed in “Demonstration of the Micro- 
wave Ice Protection Concept”, a report prepared for 
the Applied Technology Laboratory, U.S. Army Re- 
search and Technology Laboratories (USAAMRDL- 
TR-77-34) dated May 1978, and prepared under con- 
tract DAAJO2-76-C-052. While that research was pri- 
marily directed to an alternative to pneumatic and resis- 
tance heated air foil de-icers, it did report on resonant 
effects in a waveguide apparently caused by ice accre- 
tion. 

The ice detector described in the ’255 patent included 
a surface waveguide mounted adjacent on air foil in 
which ice accretion could be expected. Microwave 
energy was coupled to the surface waveguide via a dual 
directional coupler, and energy reflected back from the 
waveguide was monitored. Ice accretion, across the 
surface of the waveguide, changed the physical reflec- 
tion coefficient of the waveguide, and this could be 
detected by comparing the energy input to the wave- 
guide and the energy reflected therefrom, in the dual 
directional coupler. In the course of that research I 
noted that ice accretion apparently changed the reso- 
nant frequency of the waveguide, and ice detectors 
postulated in the cited report were based on this princi- 
ple. 

Since that time I have discovered that a device can be 
constucted for detection of icing, measuring ice thick- 
ness and detecting the rate of ice growth, by a device 
which is different from that described in either the ’255 
patent or the cited report. 

More particularly, in a first embodiment of the inven- 
tion the region to be monitored for ice accretion is 
provided with a resonant surface waveguide (or trough 
line). The waveguide is mounted flush, below the sur- 
face being monitored. Since it is non-intrusive it does 
not disturb air flow and thus may be located adjacent 
any air foil surface. Microwave energy is coupled to the 
waveguide. The microwave energy is coupled via a 
feed point. A controlled oscillator provides the micro- 
wave energy at a controllable frequency. A detector is 
coupled to the surface waveguide at a point displaced 
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from the feed point and is responsive to electrical en- 
ergy incident on the detector. A measurement device is 
coupled to the controlled oscillator for indicating the 
frequency deviation of the controlled oscillator from a 
5 quiescent frequency, and finally control means are pro- 
vided to control the frequency of oscillation of the 
controlled oscillator. In a first or open loop embodi- 
ment of the invention, an operator can operate a control 
means to vary the frequency of controlled oscillations 
10 and monitor the output of the detector, which may for 
example comprise a meter. The quiescent frequency is 
selected as the waveguide resonant frequency in the 
absence of ice. As ice accumulates the waveguide reso- 
nant frequency changes. The actual resonant frequency 
15 of the waveguide/ice combination is determined by 
adjusting the frequency of the controllable oscillator 
until the meter reads a maximum. Deviations of the 
controlled oscillator frequency, from the quiescent fre- 
quency, indicate the presence of ice, and the quantum of 
20 this deviation indicates the thickness of the ice. By cali- 
brating the measurement device, ice thickness can be 
determined directly from the meausrement device. 

If desired the measurement device can be augmented 
25 by the addition of a rate determining apparatus to indi- 
cate the rate of change of the frequency at which the 
detecting device indicates a peak. The output of the rate 
indicating circuit will therefore indicate rate of ice 
growth. 

2 q Accordingly, in an open loop embodiment the inven- 
tion provides: 

and ice detection and thickness measurement system 
comprising: 

a controlled oscillator, oscillating at a controllable 
35 microwave frequency, 

a resonant surface waveguide coupled to said con- 
trolled oscillator at a feed point, 
means to mount said waveguide in a location subject 
to ice accretion, 

detector means coupled to said surface waveguide at 
a point displaced from said feed point responsive to 
microwave electrical energy incident on said de- 
tector means, said detector means providing an 
output indicating microwave electrical energy inci- 
45 dent thereon, 

measurement means coupled to said controlled oscil- 
lator for indicating frequency deviations of said 
controlled oscillator from a predetermined fre- 
quency, and 

50 control means to control the frequency of oscillations 
of said controlled oscillator, 
whereby ice thickness is detected by noting said fre- 
quency deviation indicated by said measurement 
means when said control means controls said fre- 
55 quency of oscillation to produce a maximum de- 

tected energy at said detector means. 

In a second, or closed loop embodiment of the inven- 
tion, a feedback device is arranged to couple the detec- 
tor means to the control means. The feedback device is 
60 arranged to note for example a drop in detected energy 
by the detector means caused as a result of ice growth. 
This change in the output of the detector means pro- 
duces a signal from the feedback means to drive the 
control means in such a fashion as to change the fre- 
65 quency of oscillations of said controlled oscillator in a 
search mode to relocate that frequency which produces 
a maximum output of the detector. When the controlled 
oscillator has its frequency of oscillation changed to 
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again produce a peak out of the detector means, the 
feedback means terminates driving of the controlled 
oscillator. The frequency deviation through which the 
controlled oscillator is moved from a quiescent fre- 
quency is indicative of the ice thickness. Similar to the 5 
first or open loop embodiment of the invention, the 
measurement means can be augmented with a rate de- 
termining circuit (and including a meter to indicate the 
output of the rate determining circuit) so that the aug- 
mented measuring means indicates not only ice thick- 10 
ness, but rate of ice growth as well. 

Accordingly, in a second or closed loop embodiment 
the invention includes the apparatus recited above 
wherein the control means includes: 

feedback means responsive to said detector means for 15 
varying said frequency of oscillation to automati- 
cally determine a frequency of oscillation produc- 
ing a maximum output in said detector means. 

In contrast to the prior art cited above, and known to 
the applicant, the surface waveguide used in the inven- 20 
tion is employed in a transmission mode in that energy 
is coupled from a feedpoint to a detector, which is 
spaced therefrom, rather than in the reflection mode as 
expressed in the cited prior art. My experiments indi- 
cates that the s/n ratio of the reflection mode is so poor 25 
as to render any measurement almost useless whereas in 
the transmission mode the s/n ratio is good to excellent. 

In a preferred embodiment of the invention, one 
which for example can be employed in an operating 
aircraft such as a helicopter, to monitor ice thickness 30 
and rate of ice growth on the helicopter rotor blades, 
the feedback means is implemented in the form of a 
digital processor, for example one of the commonly 
available microprocessors now available off the shelf. 
The digital processor is provided with a program to 35 
enable it to sense detector output, via an analog-to-digi- 
tal converter, on sensing a deviation of that detector 
output indicative of a change in ice thickness, to output 
a signal to a controlled oscillator to vary the frequency 
of oscillation thereof. By controlling the output of the 40 
feedback means as a function of the output of the detec- 
tor means (using one of several conventionally known 
search strategies) the digital processor rapidly achieves 
a condition in which the detector output is maximized 
(to any given precision desired). The measurement 45 
means referred to above may take the form of analog or 
digital displays driven by the digital processor to indi- 
cate to an observer both ice thickness and rate of ice 
growth. Because of the speed with which the processor 
operates, the observer perceives, essentially in real time, 50 
information on the displays. 

Accordingly, the preferred embodiment includes two 
major components, the transducer comprising con- 
trolled oscillator, resonant surface waveguide and de- 
tector and the digital subsystem comprising processor 55 
with display output, an A/D converter coupling detec- 
tor to processor and a D/A converter coupling proces- 
sor to controlled oscillator. The transducer can be min- 
iaturized using IC’s to exhibit negligible mass. This and 
the fact that the waveguide is non-intrusive allows the 60 
transducer to be located near an air foil surface, even for 
example, a helicopter blade or engine inlet. If desired 
the processor can provide a print out of ice thickness 
and accretion rate vs. time. The processor can also store 
ice thickness and accretion rate in memory for later 65 
statistical analysis in ice or cloud studies. 

In the embodiments of the invention described above, 
only a single waveguide or transducer has been referred 
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to. In some applications it may be necessary and/or 
desirable to monitor ice at a plurality of locations within 
a single aircraft, vehicle or other structure. These needs 
are met by the use of several waveguides or transducers, 
each located according to the need and/or desire for 
monitoring ice at a variety of locations. The remaining 
apparatus, i.e. the controlled oscillator, detector means, 
measurement means and control means can be dupli- 
cated for each different waveguide or, can be time 
shared over a plurality of waveguides. Similarly, the 
second or closed loop embodiment of the invention can 
be modified by providing plural waveguides and time 
sharing some or all of the controlled oscillator, detector 
means, measurement means, control means and feed- 
back means. In a preferred form of the multi-waveguide 
embodiment of the invention, however, each of the 
waveguides is part of a dedicated transducer with a 
single processor (implementing control and feedback 
means) time shared among all transducers. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will now be further described 
in the following portions of the specification when 
taken in conjunction with the attached drawings in 
which like reference characters identify identical appa- 
ratus and in which: 

FIG. 1 is a functional block diagram of an open loop 
embodiment of the invention; 

FIG. 2 is a functional block diagram of a closed loop 
embodiment of the invention; 

FIG. 3 is an isometric illustrating a transducer, in- 
cluding waveguide, located near the surface of an air 
foil; 

FIGS. 4A and 4B are a section through the air foil at 
the location of the waveguide and a top view of the 
waveguide; 

FIG. 5 is a block diagram of an embodiment of the 
invention which has been constructed; 

FIGS. 6A-6C is a flow chart of a program for the 
processor used in the invention; 

FIG. 6D illustrates the relation between frequency 
index J and detector output; 

FIGS. 7A and 7B are plots of detector output vs. 
frequency showing changes in resonant frequency due 
to ice accretion; 

FIG. 7C is a plot of ice thickness vs. time produced 
by the inventive apparatus; and 

FIG. 8 is a plot of resonant guide wavelength vs. 
frequency in an embodiment of the invention actually 
constructed. 

DETAILED DESCRIPTION OF PREFERRED 
EMBODIMENTS 

FIG. 1 is a functional block diagram of an open loop 
embodiment of the invention which is employed for ice 
detection and measurement in connection with an air 
foil 10. As shown in FIG. 1, the air foil 10 includes a 
surface waveguide 15, which when inserted into the 
mating slot 16 in the air foil 10 exhibits a streamline 
structure. The surface waveguide 15 is shown in more 
detail in FIG. 3, at this point it suffices to note that the 
waveguide 15 includes a feed point (FP) which is sup- 
plied with microwave energy by a controlled oscillator 
17. The controlled oscillator 17 is associated with a 
control means 18 for varying the frequency of oscilla- 
tion thereof. The controlled oscillator 17 is arranged to 
oscillate at a predetermined, fixed or quiescent fre- 
quency Fq and at other times at a frequency of oscilla- 
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tion which is different from Fo as controlled by the the illustration of the air foil 10. As shown in FIG. 2, the 
manipulated shaft 18. The frequency of oscillation of waveguide 15 is driven from a VTO 29 (which replaces 
the controlled oscillator 17 is detected by a frequency the controlled oscillator 17) arid is coupled to a detector 
detector 19 which has an input coupled to the output of 21 at a point spaced from the feed point FP. The output 
the controlled oscillator 17, and an output coupled to a 5 of the detector 21, amplified in amplifier 21a is coupled 
conventional meter or other frequency indicating de- to an analog-to-digital converter (ADC) 23. The output 
vice 20. In a preferred embodiment, the meter 20 is of the ADC 23 is provided as an input to a processor 24 
arranged to indicate deviations in the frequency of the which is also coupled to associated memory 25 and I/O 
controlled oscillator 17 from the frequency Fo. As will devices 26. Processor 24 provides an output to a digital- 
be more clearly explained in subsequent portions of this 10 to-analog converter (DAC) 27 which in turn provides a 
specification, the meter 20 may be calibrated in terms of voltage to a VT driver 28. The VT driver 28 (or ampli- 

ice thickness. fier) provides an input voltage to the VTO 29 which 

Also associated with the waveguide 15 is a detector thus produces an output frequency indicative of this 

21. As is indicated in FIG. 1, the detector 21 is coupled input voltage. 

to the waveguide 15 at a point separated from the feed 15 A comparison of FIGS. 1 and 2 indicates that a feed- 
point FP, and thus the waveguide 15 operates in a trans- back device (the processor 24 and its ancillary equip- 
mission mode in that energy detected by the detector 21 ment) has been coupled between the detector 21 and the 
is energy which is transmitted through the waveguide controlled oscillator (in the form of the VTO 29). In 
15. The waveguide 15 is shorted, such as by the use of addition, the meter 20 has now been replaced by I/O 
the short S, also indicated in FIG. 1. The waveguide 15 20 device 26 which can include for example conventional 
is arranged to be resonant at the frequency Fo when digital panel meters or the like. Because of the feedback 
mated with the air foil 10. To present a streamlined provided by the processor 24, the need for the meter 22 
appearance the waveguides 15 includes an air foil sur- has been obviated, although if desired the I/O device 26 
face to mate with the air foil 10. The detector 21 is may include a read out corresponding to the meter 22. 
coupled to another conventional display such as the 25 While the processor 24 can be provided with a host of 
meter 22. The meter 22 is deflected in relation to the functions in the form of software or firmware via the 
energy detected by the detector 21. memory 25 (which may include ROM) the basic func- 

In normal operation, with the waveguide 15 installed tion of the processor 24 is to interpret variations in the 
in the air foil 10, the controlled oscillator produces an output of ADC 23 (indicative of alterations in the out- 
output at the quiescent frequency of oscillation Fo 30 put of the detector 21) and drive the DAC 27 in such a 
which is coupled to the waveguide 15. The energy is fashion as to produce a maximum or near-maximum 
detected by the detector 21 and indicated by the meter output of the detector 21. The frequency deviation of 

22. In the event that ice forms on the surface of the the VTO 29 which is required to bring the detector 21 

surface waveguide 15, the electrical characteristics of to a maximum or near-maximum condition is indicative 
the waveguide changes, and more particularly its reso- 35 of the ice accretion. When this condition is reached, the 
nant frequency changes. As a result, in the absence of a processor 24 can then interpret the digital output pro- 
change in the frequency of oscillation of the controlled ducing this condition as indicative of ice accretion or 
oscillator 17, the waveguide 15 no longer receives en- growth. Likewise, the processor 24 can also derive a 
ergy at its resonant frequency. In this off resonant con- time rate of change of ice thickness and output this 
dition, the energy intercepted by the detector 21 de- 40 (through another one of the I/O devices 26) indication, 
creases, and the meter 22 reflects this condition. The memory 25, in addition to providing the software 

An observer, observing the drop in the reading of the or firmware for operation of the processor 24 can also 
meter 22 can manipulate the shaft 18 to lower the fre- provide for temporary storage of intermediate results in 
quency of the controlled oscillator 17. Lowering the RAM. While the use of digital panel meters or the like 
frequency of the controlled oscillator increases the 45 has been indicated, the I/O devices 26 may also include 
wave length of the energy in the waveguide 15, to com- other conventional output devices such as printers, 
pensate for the change in electrical characteristics CRT’s and the like. 

caused by the ice accretion on the surface of the wave- Of course, since none of the apparatus shown in ei- 
guide 15. As the frequency of the controlled oscillator ther FIG. 1 or FIG. 2 actually measures ice thickness 
17 approaches the new resonant frequency of the wave- 50 per se, accuracy of the device depends Upon the ability 
guide 15, the energy incident on the detector 21 in- to predict the relation between frequency deviations of 
creases and the meter 22 reflects this fact. As the fre- the . controlled oscillator 17 or VTO 29 in relation to 
quency oscillation passes beyond the resonant fre- actual ice thickness. This in turn is determined by the 
quency of the waveguide 15, the energy intercepted by dimensions of the waveguide 15, the material of which 
the detector 21 again falls, reflected by movement of the 55 it is constructed, the relative locations of the feedpoint 
needle on the meter 22. The observer can then, by re- and the detector coupling, and the quiescent frequency 
versing the rotation of the shaft 18, bring the meter 22 to (Fo) of the controlled oscillator. In addition, the speed 
a maximum. During this manipulation of the shaft 18, with which the embodiment of FIG. 2 can accurately 
the frequency detector 19 is continually detecting fre- locate a new resonant frequency is determined in part 
quency deviations from the frequency Fo. Once the 60 by the software or firmware which directs the proces- 
meter 22 has again been placed in a maximum condition, sor’s search for this resonant frequency. Preferred em- 
the observer can note the condition of the meter 20 bodiments illustrating the factors which effect the fore- 
indicating the frequency deviation (Fo— F), preferably going parameters will now be described, 
with the meter 20 calibrated for the dimensions of the FIG. 3 is an isometric view of a section of a typical air 
waveguide 15 and the frequency Fo, the ice thickness 65 foil 10 within which a transducer is installed. As shown 
can be read directly off the meter 20. in FIG. 3, the waveguide 15 has a surface which gener- 

FIG. 2 is a block diagram of a closed loop embodi- ally conforms to the contours of the air foil 10, and is 
ment of the invention; the illustration of FIG. 2 omits actually mounted beneath the surface of the air foil 10. 
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The waveguide 15 is responsive to ice accumulations 
within the region defined by the uppermost surface 32 
of the waveguide 15. The surface can be protected by a 
polyurethane or polyethylene erosion shield, having a 
dielectric constant similar to the dielectric constant of 5 
the surface waveguide itself. Also shown in FIG. 3 is 
the housing 30 for the remainder of the components of 
the transducer, i.e. the oscillator 29, detector 21 and the 
associated amplifiers 28 and 21a. The equipment in 
housing 30 is coupled to the processor over a multi-con- 10 
ductor cable 31. FIG. 3 also indicates that waveguide 15 
is coupled to the housing 30 by a pair of coaxial conduc- 
tors, i.e. coaxial conductors 37 and 38. One of these 
conductors couples the VT 29 to the waveguide 15, and 
the other couples the waveguide 15 to the detector 21. 15 
The housing 30 and the components it contains can be 
miniaturized considerably by the use of microwave 
integrated circuits (MIC). 

FIG. 4A is a cross-section of the waveguide 15 and 
FIG. 4B is a top view. As shown in FIG. 4A, 35 repre- 20 
sents the surface of an air foil 10, and 33 indicates ice 
accretion above the surface of the air foil 10. The wave- 
guide 15 itself is composed of a stable dielectric material 
such as polyethylene which has approximately the same 
dielectric constant as ice. In the ice-free condition, the 25 
surface waveguide 15 is made resonant by assuring that 
its length is equal to a multiple number of half wave- 
lengths of the quiescent frequency of VTO 29. Thus, as 
shown in FIG. 4A, the length of the waveguide is an 
integral number of half wavelengths. FIG. 4A also 30 
indicates the location of the coupling 34, coupling the 
VTO 29 to the wavelength 15; more particularly it is i 
wavelength from an end wall. FIG. 4A also shows the 
detector probe 34a which is coupled to the waveguide 
J wavelength from the opposite end wall. Since this 35 
dimension is independent of ice thickness and frequency 
(guide wavelength in a resonant surface waveguide is 
constant) there is no variation in coupling due to this 
location. I have found that very light coupling is desir- 
able so a high loaded Q with very little radiation is 40 
realized. The probe is coupled to the electric field in the 
surface waveguide and therefore the coupling effi- 
ciency is dependent on the relative strength of the elec- 
tric field at the location of the probe. The coupling is 
obviously zero at frequencies below cutoff of the sur- 45 
face waveguide and starts to rise as the operating fre- 
quency is increased above cutoff. Trying to operate the 
transducer at frequencies only slightly above cutoff 
while yielding very high theoretical sensitivity will be 
impractical because the coupling will be too low. The 50 
operating frequency must be sufficiently above cutoff to 
provide adequate, albeit light coupling. The coupling is 
also a function of the height of the probe relative to the 
thickness of the dielectric; i.e. h/d, as seen in FIG. 4A. 

In the absence of ice, d=d«,. Because of the relation 55 
between coupling and the ratio h/d, the coupling de- 
creases as ice thickness increases. 

In an embodiment which has been constructed, VTO 
29 tuned in the range 5.9 to 6.57 gigahertz; the length L 
was 1.393 inches, the depth (doo) was 0.787 inches and 60 
the width W was 0.25 inches. The surface waveguide of 
these dimensions has a resonant frequency of 6.28 giga- 
hertz. 

FIG. 7A illustrates the change in resonant frequency 
of the waveguide due to increasing layers of ice. For 65 
example curve 1 indicates a resonant frequency (fi,„ at a 
frequency index fix,) in the absence of ice. As ice ac- 
cretes, the resonant frequency changes to the conditions 
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shown in curve 2 (the new resonant frequency fi, at a 
frequency index i 0 ). The change in frequency or fre- 
quency index df or di is a measure of the ice accretion. 
Curve 3 indicates a still further change in resonant fre- 
quency due to further ice accretion. 

Theoretically we could write /= \/k In (1 —df/df a ), 
where: 

df= shift in resonant frequency with ice; 
dfi,= asymptote of shift in resonant frequency with 
ice; 

k= constant; and 
I=ice thickness. 

However, I have found it more practical to use the 
relation I=l/k' In ( 1 — di/di a ), where: 
di = shift in frequency index with ice; 
di fl = asymptote of shift in resonant frequency with 
ice; and 
k'= constant. 

In the second equation, frequency index (i) replaces 
frequency (f) where frequency index is a digital count, a 
different digital count corresponding to each different 
frequency so that a measurement of the change in fre- 
quency index (di) is equivalent to a measure to change 
in frequency (df). In an embodiment which has been 
constructed, the processor generated counts in a range 
(1 to 2047) corresponding to the frequency range of 
interest. Each count was representative of a different 
frequency which could be generated by VTO 29. Thus, 
when a maximum response from the waveguide was 
determined the count (or frequency index) could be 
directly related to ice thickness. It should be apparent 
that this eliminates the necessity to actually measure 
frequency of VTO 29. It should be understood, how- 
ever, that it is within the scope of the invention to first 
vary the frequency of VTO 29 to determine the maxi- 
mum response of the waveguide (the resonant fre- 
quency) and then measure the frequency to determine 
ice thickness. The calibration process is one of deter- 
mining the values of K' and di a . 

The properties of a polyethylene surface waveguide 
may be extablished from a plot of guide wavelength 
versus frequency for various thicknesses of the surface 
waveguide as illustrated, in FIG. 8. FIG. 8 shows this 
relation between 1 and 2 Ghz for the TE; mode. This is 
obtained by conventional techniques on a trough line 
whose dimensions are given above. Similar plots made 
over other octave bands such as 2 to 4 and 4 to 8 Ghz, 
etc., reveal that all plots are identical if the scales are 
related by a normalization factor p so that FIG. 8 may 
be used to represent any octave band given the normal- 
ization factor p. 

Referring to FIG. 8, the surface waveguide of thick- 
ness d is made resonant by making its length equal to an 
integral number of half wavelengths. Thus, from FIG. 8 
when dp =8 and f/p= 1.66 Ghz (point “a” FIG. 8), the 
surface waveguide is resonant when its length 
p\g=14.2 cm. This is defined here as the quiescent 
resonant frequency, fix,. Of course, guides of different 
thicknesses could use different fix,, and while polyethyl- 
ene is a preferred guide material, other materials could 
be used with corresponding changes in fix,. 

The growth of ice layers, assuming the ice has ap- 
proximately the same dielectric constant as the base 
dielectric, has the effect of increasing the thickness d of 
the surface waveguide. As seen from FIG. 8, this has 
the effect of decreasing the wavelength for a given 
frequency and destroying the resonance. The condition 
of resonance is restored by increasing the wavelength 
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by decreasing the frequeny until the wavelength is one 
again equal to the guide wavelength that existed in the 
ice-free condition. The turning process is illustrated in 
FIG. 8 as the locus of points along the horizontal line of 
constant X go. 

While a polyethylene waveguide has been used, this 
technique allows waveguides of different materials to 
be used. Preferably any other waveguide should have 
dielectric constants near to that of ice. 

In testing the transducer whose dimensions are given 
above, changes in ice thickness was simulated by adding 
layers of polyethylene tape (since the dielectric constant 
of ice and polyethylene are close to each other). A pen 
recording, FIG. 7B, was produced to show the shift in 
resonant frequency for tape thickness of, respectively, 
0, 12, 24, 36, 48, 72, 84, 96 and 108 mils. At zero thick- 
ness the resonant frequency was 6.29 gigahertz each 
corresponding change in frequency was 0.05, 0.95, 
0.120, 0.142, 0.185, 0.250, 0.227 and 0.253 gigahertz. 

The curve in FIG. 7C illustrates a real time pen re- 
cording of output indications of ice thickness as a func- 
tion of time produced in accordance with the invention. 

FIG. 5 is a block diagram of the apparatus which had 
been constructed to produce the curves shown in FIG. 
7C. Shown in FIG. 5 is VTO driver 28 which drives the 
VTO 29. In the embodiment actually constructed VTO 
29 was the Avantek Model 8580 tuning from 5.8 to 6.6 
gigahertz. Detector 21 is a crystal detector, its output 
after amplification at 21a, is input to a two-channel 
analog-to-digital converter 23; one channel with a gain 
of one and the other channel with a gain of three. The 
converted signal output of the transducer in both chan- 
nels are made available to the microprocessor 24; in the 
embodiment constructed this was a Texas Instrument 
Tm990-101M using a TMS9900 16-bit microprocessor 
and 16K-bytes RAM. Of the apparatus shown in FIG. 
5, the printer 26a, keyboard 266 and the digital-to- 
analog converter, scope and plotter 26 d, while useful 
for development, are not necessary to system opertion. 
On the other hand, the I/O device 26c representing a 
TM990/2 10 input/output expansion board provides a 
driver for two output ports driving a first panel display 
device (PDD) displaying ice thickness, and a second 
panel display device (PDD) displaying rate of ice 
growth. 

In order to operate automatically, once calibrated, 
the microprocessor drives the D/A output with a par- 
ticular number representing a frequency index. This is 
converted first to analog form (by the digital-to-analog 
converter 27 and the amplifier 28) and then as a voltage 
is converted to a frequency by the VTO 29. As a result 
of this particular frequency, the detector 21 detects a 
particular voltage which is amplified (by the amplifier 
21a) and then converted to digital form by the two- 
channel analog-to-digital converter 23. Depending 
upon amplitude, one of the two digital signals is selected 
by the microprocessor. The microprocessor reads this 
quantity as the response of the transducer to this partic- 
ular frequency. 

The computer then emits another frequency index 
number and processes the return signal looking for a 
maximum return signal. When a maximum is found, the 
index which produced it is the resonant frequency in- 
dex. 

The program which is disclosed below uses ML for 
resonant frequency index (i 0 ). The quiescent frequency 
index (i m ) is referred to as CA. Thus, the difference 
CA —ML =DX is the change in frequency index (di) 
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caused by ice accretion. A function of this difference 
DC is then related to ice thickness. . 

In any particular scan, when ML is found, it is used 
by the computer to determine the range of indices J to 
5 be presented next, J1 to J2, and the process of searching 
for the peak between J1 and J2 is repeated. In the em- 
bodiment actually constructed, the sampling rate is six 
thousand frequencies per second. Every one-half sec- 
ond for the first thirty seconds and every second there- 
10 after, at the completion of a sweep, the latest values 
DX, the time, the ice thickness data, and the amplitude 
of the peak MX, are stored in an array in computer 
memory. At the same time, the ice thickness data and 
accretion rate are converted to decimal and presented 
15 to the two panel digital displays. 

FIG. 6D is a simplified curve of a single sweep plot- 
ting Aj (the returned voltage from the detector 21) and 
the frequency index J. In this processing, J is decre- 
mented from J1 to J2. MX is the highest value of Aj 
20 found. MJ is the value of J at which MX is first encoun- 
tered. MK is the value of J at which MX is last found 
(which may be equal to MJ). ML is the average 
(MJ+MK)/2. After a sweep, 71= ML +48 and 
72 =ML— 48 unless: 

25 MX— 72 <3, then 71=72+16; 72=71-96, or 
71— ML<3, then 72=71-16; 71=72 + 96, or 
ML > CA, then 71 = CA + 8; 72 = CA — 88, or 
MX <500 or OS <8, then 71=G4; 72=1. 

If MX < 500 and OS<8, then OS is incremented by 8. 
30 If MXy 2000 and OS>7, then OS is decremented by 8. 
OS is the amplification flag, if not less than 8 the high 
gain channel is used, Aj is amplified by 3, otherwise a 
unity gain channel is used. 

A flow chart for the program used in the embodiment 
35 we actually constructed is shown in FIGS. 6A-6C. 
Referring first to FIG. 6A, functions F2 and F4 ini- 
tialize the processor including loading the program, 
initializing registers, storing the assembly program, and 
clearing arrays. Function F5 runs a test sweep through- 
40 out the entire range (in our embodiment frequencies 
indices 1 to 2047 corresponding to the entire frequency 
range of the VTO). The parameters displayed are the 
highest return amplitude from the detector 21 (MX) and 
the average of the indices over which that amplitude 
45 was returned (ML). Function F6 allows an operator 
input to initiate calibration. Assuming that calibration is 
desired, function F7 is performed which sweeps 
through a smaller frequency range (in our embodiment 
64 indices centered about ML). The sweep is performed 
50 thirty times and a parameter ML is the frequency index 
at which the highest returned amplitude is found. A 
parameter Sm is the sum of the parameter ML over all 
thirty sweeps. 

Function F8 then computes a number of parameters, 
55 CA is the average ML over the thirty sweeps and is 
considered the ice-free quiescent resonant frequency 
index. An ice trigger parameter IT is also computed, the 
parameters MX, CA are stored into the arrays and the 
parameters CA and IT are displayed on the panel dis- 
60 play meters. Function F9, in our embodiment is then 
performed; however, in the absence of a CRT display, 
function F9 is not required. Function F10 allows the 
operator to request recalibration. Assuming recalibra- 
tion is not effected, a loop of functions F9 to Fll is 
65 performed as long as the difference between CA and 
ML is less than the ice trigger IT. When that inequality 
is no longer satisfied, icing has been detected and pro- 
cessing skips thrdugh point B (FIG. 6B). Function F12 



4 , 470,123 


11 

sounds an alarm and starts a timing clock. Function F13 
sets up new parameters J1 and 32 (the limit within 
which the frequency index is swept) as a count of 96 
centered on ML. Function F15 performs a single sweep 
and returns with four parameters, MX is the highest 5 
output amplitude of the detector 21 within the sweep, 
MK is the index value at which MX was last found, MJ 
is the index value at which MX is first detected, and ML 
is the average. This is effected by a subroutine PEAK, 
shown in more detail in FIG. 6C. Thereafter, function io 
F16 is effected. The adjustments to Jl, J2 and OS, based 
on comparisons between JK, J2; Jl, ML; ML, CA; MX, 
500; OS, 8; MX, 2000 and OS, 7 are as described above. 

Functions F15, F16 are repeatedly performed until 
the clock, checked at F28 is expired. When the clock 15 
expires, function F30 computes the difference between 
the frequency index (ML) and the calibration index 
(CA). If DX is less than zero (determined at function 
F31), then function F32 is performed to change the 
sweep parameters. TH (ice thickness) is set to zero. On 20 
the other hand, if DX is greater than zero, then function 
F33 computes TH from DX. 

The ice thickness TH is computed from DX based on 
the fifth order polynomial 2C/DX' where the coeffici- 
ents are: 25 

0.426986341 
-0.00009568097 
3.50247865 XlO- 7 
-3.49928882x10-1° 

1.25803687X10-13 
respectively. 

Thereafter, function F34 determines the time at 
which measurement was effected and stores the param- 
eters MX, DX and TH in the array. Thus, as functions 
F15-F33 are repeatedly performed the array created by 
F34 retains more and more data. Smoothed ice thick- 
ness TT and rate TK are obtained from the six most- 
recent readings of TH and fitted to a least squares 
straight line. Thus, F39 checks to see if six or more 
entries have been accumulated. If so, F40 is checked to 
see if an initial fit has been obtained. If not, F42 is per- 
formed. On the other hand, if curve fitting has been 
effected, the F43 merely updates sums based on new 
data. In either event, F46 converts binary to decimal 
and displays. Thereafter, F41 again checks the timer to 
see if a reduced sampling rate can be effected by com- 
paring elapsed time to a fixed quantity (e.g. 30 seconds). 
Depending on the outcome of the comparison, F45 may 
be effected. Thereafter, F44 checks to see if the array is 
full. If not, the loop is reentered. If so, printing of data 
is implemented. 

In an installation which did not include a printer, 
once icing was detected, the loop of functions F15-F46 
is continually performed to update the display of ice 
thickness and rate of growth of ice thickness. 

The subroutine PEAK, referred to at function F15 
(FIG. 6B) is shown in more detail in FIG. 6C. The 
subroutine is entered with the parameters Jl and J2. 
These parameters are determined, in the first pass at 
function F13, and in succeeding passes, at function F16. 

In any event, F15.1 moves the parameter Jl to the regis- 
ter J. Function F15.3 loads a suitably low quantity into 
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tion of the analog-to-digital conversion, either function 
F15.12 or function F15.14 stores the quantity A. There- 
after, function F15.15 decrements J and function F15.16 
determines if J is less than J2. If it is, the subroutine 
loops back to function F15.4 and outputs the new pa- 
rameter J and awaits the return signal. 

On the other hand, the returned amplitude would 
normally be greater than MX, since the initial quantity 
is a low quantity. Thereafter, function F15.7 updates 
MX, and function F15.8 updates the register MJ with 
the present parameter J. Function F15.9 updates the 
register MK with the present parameter J and the OS 
flag is checked again at function F15.10. 

This processing is continued until a parameter A is 
detected which is equal to MX. At that point, function 
F15.6 skips the processing to function F15.9. Skipping 
function F15.7 and F15.8 maintains previous quantity 
MX and previous quantity MJ. This is significant since 
these values will normally not be changed so long as the 
quantity A is either maintained level or decreases. In 
this fashion, the subroutine identified MJ as the first J 
parameter at which the largest return amplitude was 
detected, and MX as that amplitude. 

This processing continues until function F15.5 deter- 
mines that a return quantity A is less than MX. As a 
result, in addition to skipping functions F15.7 and F15.8, 
function F15.9 is also skipped, thereby retaining MK as 
the last parameter J at which the highest return value 
MX was detected. 

This processing then continues until functions F15.17 
determines that J is not less than J2. Thereafter, func- 
tion F15.18 computes the average ML and that process- 
ing terminates and returns to function F16. 

Reference again to FIG. 6D illustrates the expected 
relation between the quantitites Jl, J2, MJ, MK and 
ML. 

A review of FIGS. 6B and 6C indicates that function 
F15 and particularly function F15.4 outputs a sequence 
of frequency indices to sweep the VTO 29 in frequency 
within a predetermined range (Jl to J2). The functions 
F15.5 through F15.9 assist in locating the returned peak 
(i.e. MJ and MK) and isolate the peak digital signal MX. 
Function F15.18 in determining ML, in conjunction 
with function F16, in determining Jl and J2, determine 
the predetermined frequency range which will be 
swept. Finally, F33 converts the frequency index differ- 
ence DX into an ice thickness signal TH. As explained, 
TH is smoothed by functions F39-F42 to derive 
smoothed ice thickness TT and ice rate signal TK. 

Those of oridinary skill in the art will understand that 
the invention may be used in conjunction with a con- 
ventional de-icer, which for example can be turned on 
or initiated once the invention indicates substantial ice 
accretion. After a period of operation of the de-icer, of 
course, the ice will be removed. In that event, the func- 
tions F31, for example, will indicate lack of ice. At that 
point, the program can be returned to the loop of func- 
tions F10 and Fll to await further icing, in a manner 
that should be apparent to those skilled in the art. 

We claim: 


a register MX. Function F15.4 outputs the parameter J l. An ice detector and thickness measurement system 
(at this point equal to Jl) and starts the analog-to-digital comprising: 
conversion. The resulting quantity is then compared 65 a controlled oscillator, 

with the quantity in MX. First, at function F15.5, if A is a shorted surface waveguide coupled to said con- 

less than MX, processing skips to function F15.10 trolled oscillator at a feed point with a short lo- 

where the amplification flag is checked. Depending on cated a fixed distance from said feed point, 

the condition of the amplification flag, and on comple- 
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means to mount said waveguide in a location subject 
to ice accretion, 

detector means coupled to said surface wave guide at 
a point displaced from said feed point and respon- 5 
sive to electrical energy incident at said detector 
means for indicating electrical energy incident 
thereon, 

measurement means coupled to said controlled oscil- 
lator for indicating frequency deviation of said 
controlled oscillator from a predetermined fre- 
quency, and 

control means to control a frequency of oscillation of 
said controlled oscillator, 

whereby ice thickness is detected by noting said fre- 
quency deviation of said measurement means when 
said control means controls said frequency of oscil- 
lation to produce a maximum detected energy at 
said detector means. 

2. The apparatus of claim 1 in which said control 
means includes: 

feedback means responsive to said detector means for 
varying said frequency of oscillation to automati- 
cally determine a frequency of oscillation produc- 
ing a maximum of detected energy by said detector 
means. 

3. The apparatus of claim 2 in which said measure- 
ment means, said control means and said feedback 
means comprise a microprocessor with an output port 
and an input port, 

a digital-to-analog converter coupled between said 
output port and said controlled oscillator, 

an analog-to-digital converter coupled between said 
detector means and said input port, and wherein 
said microprocessor includes: 
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frequency sweeping means to provide at said output 
port, digital signals to sweep said controlled oscil- 
lator through a predetermined frequency range, 
peak following means responsive to signals at said 
input port to isolate a peak digital signal and to 
determine said predetermined frequency range, 
and 

conversion means to develop an ice thickness signal 
from said peak digital signal. 

10 4. The apparatus to claim 1 wherein: 

said controlled oscillator is a voltage-tuned oscillator 
and wherein said measurement means and said 
control means comprise a microprocessor with 
output and input ports, 

15 an analog-to-digital converter coupled between said 
detector means and said input port, 
a digital-to-analog converter and an amplifier serially 
connected between said output port and said volt- 
age-tuned oscillator, and wherein: 

20 said microprocessor includes: 

means responsive to signals at said input port to iso- 
late a peak digital signal, and 
conversion means for converting said peak digital 
signal to an ice thickness signal. 

25 5. The apparatus of claim 1 wherein: 

said shorted surface waveguide is a dielectric trough 
line with a dielectric constant near to a dielectric 
constant of ice. 

6 . The apparatus of claim 5 in which said trough line 
30 is polyethylene of a length resonant at said predeter- 
mined frequency. 

7. The apparatus of claim 6 in which said feed point is 
located one-fourth wavelength from an end of said 
trough line. 

8. The apparatus of claim 7 wherein said trough line 
has a depth h and therein said feed point is located a 
distance d below said waveguide surface wherein a 
ratio d/h provides light coupling in a TE1 mode. 
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